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Cet exposé est consacré aux avancées mathématiques récentes sur un problème célèbre 
de l'astrophysique, la stabilité de modèles galactiques. La question se formule très sim- 
plement : si l'on considère un ensemble d'un très grand nombre d'étoiles en interaction 



> 

^ ■ corps 



gravitationnelle^ 1 ^ dont la cohérence est assurée par leur attraction réciproque et que 
l'on considère en première approximation comme un système fermé, quelles sont les 
répartitions statistiques stables au cours du temps? Autrement dit, quelles sont les 
configurations de galaxies observables dans notre univers? On néglige ici les effets re- 
lativistes et on se place dans le cadre de la mécanique classique. 

Ce problème semble à première vue tout autant insoluble que le problème à N corps 
de Newton. Comment formuler des prédictions à long terme sur un système de 10 11 



(2) 



alors même que l'on ne sait pas résoudre de manière satisfaisante le problème 
à 3 corps ! C'est ici que la mécanique statistique entre en jeu. En suivant les idées de 
^ . Maxwell et de Boltzmann, on peut tenter de décrire de manière statistique l'évolution 

de nos N corps lorsque N tend vers l'infini et que les corps sont suffisamment « peu 
corrélés », à travers une équation aux dérivées partielles non-linéaire sur la répartition 
d'un corps typique. 

Cette approche a d'abord été appliquée aux cas de gaz collisionnels pour donner 
la célèbre équation de Boltzmann. Cependant, les collisions entre étoiles dans une ga- 
laxies sont quasi-absentes et l'interaction se fait essentiellement à distance via le champ 
gravitationnel. Dans les années 1930, Vlasov a découvert comment effectuer la limite 
N — y +oo dans ce cas « non-collisionnel » afin d'obtenir des équations aux dérivées par- 
tielles non-linéaires dites de « champ moyen ». La plus célèbre d'entre elles est l'équation 
de Vlasov-Poisson, dont nous allons étudier ici la version gravitationnelle dHJ)-©- Cette 
équation de transport non-linéaire décrit avec une excellente précision l'évolution de 
systèmes stellaires sur de grandes échelles de temps. 



1. Le rôle joué par les planètes dans la dynamique est négligé au premier ordre car leurs masses 
sont bien plus petites que celles des étoiles. 

2. C'est l'ordre de grandeur du nombre d'étoiles dans notre galaxie. 
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Ainsi, lorsque le nombre de corps est grand et que les corrélations sont faibles, on peut 
espérer formuler des prédictions de stabilité à partir de l'étude de l'équation de Vlasov- 
Poisson. C'est le physicien russe Antonov |An14 IAn2] qui, le premier, découvre com- 
ment résoudre le problème de la stabilité, mais uniquement dans un cadre linéarisé. Il 
démontre ainsi la stabilité linéarisée des modèles galactiques sphériques et monotones en 
l'énergie microscopique pour des petites perturbations. Cependant, l'équation de Vlasov- 
Poisson est réversible en temps, non dissipative, et ses solutions montrent des oscillations 
en temps grand ; rien ne garantit a priori que l'étude de stabilité linéarisée d'Antonov 
implique la stabilité non-linéaire recherchée. C'est ce problème mathématique que nous 
appellerons conjecture de stabilité non-linéaire à la Antonov. Durant les cinquante 
dernières années, l'analyse mathématique des équations cinétiques de Boltzmann et 
Vlasov a connu un fort développement. Cette conjecture de stabilité non-linéaire a été 
récemment démontrée par Lemou, Méhats et Raphaël [LMR61 [LMR8| lLMR9j . suivant 
des travaux précurseurs de Dolbeault, Guo, Hadzic, Lin, Rein, Sânchez, Soler, Wan, 



Wolansky ainsi que Lemou, Méhats et Raphaël [WâTl PWôH IWâ2l ICûTl IGÛ21 ICûâl [Gû4l 
IGRT1 IGR21 IDSS1 ISSl [Hal IGR31 ICLl iLMRTl ILMR31 [LMR4] . 



Dans cet exposé, nous proposerons tout d'abord dans la section 1 une introduction 
mathématique au problème, en expliquant l'origine du système de Vlasov- Poisson gravi- 
tationnelle à partir du problème à N corps. Nous rappellerons ensuite dans la section 2 
les principales propriétés mathématiques de ce système d'équations. Puis nous aborde- 
rons dans la section 3 la résolution du problème linéarisé. Enfin nous traiterons dans 
la section 4 de la question de la stabilité non-linéaire. Nous retracerons l'histoire des 
différentes méthodes mathématiques développées pour attaquer ce problème, puis nous 
détaillerons le théorème central du travail |LMR9] de Lemou, Méhats et Raphaël, en 
donnant un schéma détaillé de preuve. Nous conclurons finalement avec des commen- 
taires et questions ouvertes. 

L'auteur remercie Y. Guo, P.-E. Jabin, M. Lemou, F. Méhats, Z. Lin, P. Raphaël 
et J. Soler pour les échanges par courriel ou discussions durant la préparation de cet 
exposé, ainsi que S. Martin pour ses relectures et commentaires sur ce texte. 



1. PROBLÈME À N CORPS ET SYSTÈME DE VLASOV-POISSON 

1.1. Le problème à N corps de Newton 

On écrit le problème à N corps, pour des interactions binaires et sans champ extérieur, 
et en normalisant toutes les masses à 1 : 



Viit) 



dH 

dvi 



(X(t),V(t)) 



(1) 



1 < i < N, 
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dH 
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où tjj est le potentiel de l'interaction. Ces équations sont écrites sous forme hamiltonienne 
pour le hamiltonien 

(2) H(X,V) = J2 l -^f + Y,^ x *-^ 

i=l i<j 

où X = (xi, . . . , aîjv), V = (vi, . . . , va) et chaque Xi, Vi G M 3 . 

Même si ce système d'équations différentielles ordinaires non-linéaires couplées peut 
être aussi bien utilisé pour décrire l'évolution d'objets infiniment petits, comme des 
molécules ou infiniment grands, comme des étoiles, nous parlerons ici de point de vue 
microscopique et trajectoriel et nous appellerons i/j le potentiel d'interaction microsco- 
pique et H le hamiltonien microscopique. 



1.2. La limite de Vlasov ou limite « de champ moyen » 

Si l'on considère l'évolution de la distribution f( N ' = f^ N \t, X, V) de nos N corps 
dans l'espace des phases, des positions et vitesses, on obtient V équation de Liouville à 
N corps, qui est le point de vue statistique sur ce système : 

1=1 x ' 

où {•, •} désigne le crochet de Poisson sur (IR 3 )^ x (IR 3 )^ : 

Les solutions sont données par f^ N \t, X, V) = f( N \S- t {X, V)) où S t désigne le flot 
solution des équations de Newton. 

Cette résolution par la méthode des caractéristiques implique la conservation du signe 
de /W : 

f {N) (0, X, V) > f iN) (t, X, V) > 

mais aussi, du fait que le jacobien de S t vaut 1 pour tout t K 3 q de la masse 

/ fW(t,X,V)dXdV = [ f {N) (0,X,V)dXdV 

i(R3 x R3)iV J(R3 X R3)iV 

et plus généralement de toute < fonctionnelle de Casimir » 

/ C (f m {t,X,V)) dXdV = [ C(f w (0,X,V)) dXdV 

pour C G C 1 (M + ,IR + ), C(0) = 0. Une fonctionnelle de Casimir désigne pour un système 
hamiltonien une fonctionnelle qui annule le crochet de Poisson avec toute autre fonc- 
tionnelle le long des trajectoires (voir par exemple |Arl[ Rr2] ). 



3. Ce qui est bien sûr une version du célèbre théorème de Liouville dans ce contexte. 
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Afin de comprendre l'évolution à travers un système réduit, on introduit la distribu- 
tion d'une particule / = : 

f(t,x,v):— / f (t, X, V) dx 2 dx3 . . . dx^v dv 2 dt>3 . . . dv n 

qui est donnée par la marginale de la distribution complète à N corps f^ N \ On calcule 
alors l'équation vérifiée par cette distribution à une particule 

(4) ?L + V .vj-(N-1) ! V x ij>(x-y)-V v f®(x,y,v t vJdydv. = 0. 

Cette équation dépend bien sûr de la seconde marginale (distribution à deux parti- 
cules) 



f i2) (t,xi,x 2 ,vi,v 2 ) := / f {N) dx 3 dx A . . . dx N dv 3 dv 4 . . . dv N 

J(R3xR3)iV-2 

qui contient de l'information sur les corrélations du système. De la même façon, on 
pourrait définir les fc-marginales f^ k ' pour 1 < k < N, et l'ensemble des équations 
couplées sur toutes ces marginales constitue la hiérarchie BBGKY (Bogoliubov-Born- 
Green-Kirkwood-Yvon) . 

La limite de champ moyen consiste : 

- d'une part à considérer que chaque interaction binaire est d'ordre 0(1/N), soit 

ï/j = i/j n = —, 

ce qui signifie que l'influence d'une étoile sur une autre étoile est négligeable à 
l'échelle de l'ensemble de la galaxie, mais que l'action de l'ensemble de la galaxie 
sur une étoile donnée est d'ordre 0((iV— 1)/N) = 0(1); 

- d'autre part à considérer que 

(5) / (2) ~ / x / lorsque N -> oo, 

hypothèse que l'on appelle propriété de chaos dans ce contexte, et qui traduit des 
corrélations faibles dans le système. 
Si l'hypothèse §5§ reste vraie au cours du temps lorsqu'elle est vérifiée au temps 
initial, on parle de <c propagation du chaos ». Démontrer cette propriété est un problème 
mathématique ouvert dans le cas de l'interaction gravitationnelle (mais également pour 
l'interaction de Coulomb répulsive entre particules de même charge) ; nous renvoyons à 
|Dot IBIi] IHJlt IHJ2] pour des résultats partiels sur ce sujet. 

Sous réserve de l'hypothèse d'échelle sur tp^ et de la propagation du chaos, on obtient 
l'équation fermée suivante sur la distribution / dans la limite N — > +oo : 

(6) |+fVJ- / Vj(x - y) ■ V v f(x, v) f(y, w) dy dw = 0, 

at JR3 X R3 

que l'on peut réécrire en 

(7) ^ + v ■ VJ - H a Vj(x - y) H 3 f(y, v,) dvA dy \ ■ V v f(x, v) = 0. 
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1.3. Le système de Vlasov-Poisson gravitationnel 

Lorsque le potentiel microscopique correspond aux interactions gravitationnelles 
ip(x) = — l/(47r|x|), on obtient ainsi le système de Vlasov-Poisson gravitationnel : 

f dtf + v- V x f - V x <j) f ■ V v f = 0, te R+, x G M 3 , v G M 3 , 

(8) 

l f(t = 0,X,V) = f in (x,v) 

où l'on définit le champ moyen gravitationnel 0/ (dépendant de la fonction /) par 

(9) pf(t,x):=f f(t,x,v)dv et <j> f (t, x) :=-——* Pf . 

J R 3 4ii\x\ 

Observons que le champ <fif se résout par V équation elliptique de Poisson 

(10) A x( p f = Pf 

ce qui confère son nom au système de Vlasov-Poisson. 

On voit ici le caractère fondamental du système de Vlasov-Poisson. Dans le cas 
d'interactions à distance (par opposition aux mécanismes de collision), il est à la 
mécanique classique statistique ce que les équations de Newton sont à la mécanique 
classique. 



1.4. Le problème de stabilité d'un point de vue statistique 

On peut maintenant reformuler la question de la stabilité des galaxies dans le cadre 
du système de Vlasov-Poisson. Cela rejoint l'une des premières questions que l'on se 
pose face à une équation aux dérivées partielles : quelles sont les solutions stationnaires 
de ce système et, parmi ces dernières, lesquelles sont stables ou instables? 

De nouveaux aspects mathématiques différents de l'analyse des équations de Newton, 
et sur lesquels nous reviendrons plus loin, surgissent alors : quelles notions de solution 
utiliser pour l'équation ? Quels espaces fonctionnels utiliser ? Quelle régularité pour ces 
solutions ? 

Puisque nous sommes partis des équations de Newton, il est légitime de se de- 
mander quel serait l'impact en retour d'un théorème de stabilité sur le problème à 
N corps. C'est une question difficile et mal comprise à l'heure actuelle. Nous pouvons 
néanmoins avancer deux principes : d'une part, cela dépend fortement des réponses 
données aux problèmes ouverts de limite de champ moyen et propagation du chaos et, 
d'autre part, cela dépend vraisemblablement également de l'espace fonctionnel dans 
lequel les résultats de stabilité sont obtenus. Dans cet exposé, nous nous intéresserons 
à des voisinages de stabilité pour des espaces des Lebesgue ; ces derniers sont peut-être 
trop peu réguliers pour espérer en déduire des informations sur le problème à N corps. 
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2. PROPRIÉTÉS DU SYSTÈME DE VLASOV-POISSON 

Nous allons tout d'abord rappeler les propriétés mathématiques élémentaires des 
équations (JE])-®. 

2.1. Énergie microscopique et hamiltonien 

Les équations de Newton préservent le hamiltonien microscopique H défini par §2§ : 

VtGl, H(S t (X,V)) = H(X,V). 

En intégrant ceci contre la distribution initiale de nos N particules (équation de 
Liouville) et en divisant par N (pour éviter la divergence de l'énergie lorsque N tend 
vers l'infini), on obtient la conservation statistique du hamiltonien microscopique 

vteR , / H{X A ; V) f {N \t,x,v)dXdv=[ H{X A ; V) f {N \0,X,V)dXdV. 

Dans la limite de champ moyen, cette égalité implique par calcul élémentaire la conser- 
vation au cours du temps de la quantité 

W)= / l -^f(t,x,v)dxdv + [ ^Af( t ,x,v)dxdv 

JM 3 xR 3 1 JR 3 xR 3 * 

que l'on appellera hamiltonien (macroscopique) du système de Vlasov-Poisson. Remar- 
quons le facteur 1/2 devant 0/ qui rappelle que notre équation sur la distribution à une 
particule / est la « trace » d'un système à TV particules avec interaction binaire. Cette 
fonctionnelle se réécrit en utilisant l'équation de Poisson (TTUj) sur 0j : 

(H) «(/)=/ l -^f(t,x,v)dxdv- [ \Y^MlÉl àxàv . 

jR 3 xl 3 1 Jm 3 xr 3 1 

On voit ici la possibilité de transferts d'énergie entre énergie cinétique et gravitation- 
nelle. De plus, le caractère attractif de l'interaction se traduit par le fait que ces trans- 
ferts peuvent « diverger » : la fonctionnelle % conservée est la somme de deux termes 
de signe opposé qui peuvent potentiellement diverger tout en s'équilibrant, ce qui est 
une source de difficultés mathématiques. 

Notons enfin que l'équation de Vlasov de champ moyen ([8]) peut être interprétée 
comme l'équation de Liouville sur une particule associée à V hamiltonien microscopique 
de champ moyen 

M 2 

E = E f (t,x,v) = y+^,x) 

sous la forme 

d t f + {f,E f } = 
où {-, ■} désigne ici le crochet de Poisson sur K 3 x R 3 : 
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Remarquons qu'il n'y a plus le facteur 1/2 devant l'énergie potentielle dans la définition 
de Ef : la cause en est que l'on considère ici une particule évoluant dans un champ moyen 
donné en « oubliant » la non-linéarité. 

2.2. Fonctionnelles de Casimir et équimesurabilité 

L'équation de champ moyen (JSJ) est associée aux solutions (dites caractéristiques) du 
système différentiel 

x'{t) = v{t) 



dE f 



dv 

v'it) = -V x <j> f (x{t)) 



dE f 



dx 

Ces courbes dépendent bien sûr de la solution elle-même, mais cela montre que 
l'équation d'évolution (jSJ) est une dynamique de réarrangement de la distribution /, 
qui préserve le signe, mais aussi toutes les fonctionnelles de Casimir 

C(f(t,x,v)) dxdv = / C(f(0,x,v)) dxdv 

pour C G C 1 (M + ,]R + ), C(0) = 0. Cet ensemble de lois de conservation peut être 
synthétiquement exprimé par la propriété suivante : la solution reste pour tout temps 
équimesurable à sa donnée initiale 

Vt > 0, Vs > 0, \{\f{t,x,v)\ > s}\ = \{\f in (x,v)\ > s}\ G K+ U {+00} 

où \0\ désigne la mesure de Lebesgue d'un ensemble mesurable 0cK 3 x R 3 . 

Nous sommes donc en présence d'un système dynamique en dimension infinie 
possédant une infinité non dénombrable de contraintes : le hamiltonien macroscopique 
%(f) et l'équimesurabilité. 

2.3. Lien entre invariants et hamiltonien : une inégalité d'interpolation clef 

Comme nous l'avons vu, l'énergie totale du système TL(f) est la différence de deux 
termes positifs que nous noterons 

^(/) = Hcinif) — Hpot(f)- 

Une difficulté évidente sera alors de contrôler ces deux formes d'énergie séparément. A 
cette fin, il est possible d'interpoler l'énergie potentielle à partir de l'énergie cinétique 
et d'un espace de Lebesgue L p , dont la conservation est assurée par l'invariance des 
fonctionnelles de Casimir. On peut ainsi montrer les inégalités fonctionnelles suivantes. 

PROPOSITION 2.1. — Pour tout p G [1,+oo], on a l'estimation d'intégrabilité de la 
densité spatiale 

" P/ "l§^T(R3) - C \\f\\lJ'~(M.6 ) 'Hcin(f) 1 

avec dans le cas limite p = +00 

ll^/H rirrcai - C \\f\\l^(w6 ) 'Hcin(f)^, 



3p-3 
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pour une certaine constante C > 0. 

On en déduit le contrôle suivant de l'énergie potentielle pour p > p c = 9/7 

Hpot{f) < C|l/llLi î (M6)||/|| i p î) ( R 6)'Wcm(/) 2 

pour une certaine constante C > 0. 

Démonstration. — Pour p > 1 et R > 0, on décompose l'intégrale de la densité spatiale 
Pf en 



fdv = fdv + fdv 

J\v\<R J\v\>R 

(p-1) 

<\\f\\r*([ àv) P +R~ 2 I f\v\ 2 dv 

\J\v\<R J 

<\\f\u(^f-) v +R~ 2 



En optimisant le paramètre R, on obtient 



(3p-3) 

2£ / P \ (5p-3) 



fdv<C\\f(x,-)\\^ 3) [ / f\v\ 2 dv 



et on en déduit par inégalité de Cauchy-Schwarz 




(3p-l) 

/ „ \ ( 3 P~ 3 ) \ (5p-3) 

2 P / /" \ (3p-l) 



(3p-3) 



ce qui conclut la preuve de la première inégalité. 

On veut maintenant contrôler le champ 0/. En raisonnant heuristiquement, l'inégalité 
de Sobolev et la régularité elliptique de l'équation de Poisson impliqueraient 

||V x 0/||l2(1R3) < CHV^/ 11^1,6/6(^3) < C||p/|| L 6/5( K 3). 

La fonction V x 0/ n'est pas dans l'espace iy 1,6 / 5 (IR 3 ) ; cependant on a bien l'inégalité 
11^0/11^2^3) < (^1^/11^6/5^3) en utilisant la formule V x 0/ = (x/(47r|x| 3 )) * p/ et 
l'inégalité de Hardy-Littlewood-Sobolev \LL\ Theorem 4.3]. De manière alternative, 
il est possible de montrer que V x 0/ G 

(1 'espace de Sobolev homogène) au moyen 
du contrôle sur le laplacien de 0/ et de la théorie de Calderôn-Zygmund, puis d'utiliser 
l'inégalité de Gagliardo-Nirenberg-Sobolev pour conclure. 

En combinant ces deux précédentes inégalités, on obtient 

Upot(f) = W^x4>f\\h(M.3) < C||P/||^ 6 /5( R 3)- 
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Afin d'estimer ce membre de droite on effectue l'interpolation suivante : pour 
p > Pc = 9/7, on a (5p - 3)/(3p - 1) > 6/5 et 

(7p-9) (5p-3) 
Il || s || || (12p-12) || || (12p-12) 

l|P/llL6/5(R3) S ||P/|| L 1( R 3) HP/H (5 P -3) 

L(3p-i) 

(7p-9) (5p-3) 
< Il /Il ( 12 P- 12 ) Il n .11 (12P-12) 



L 1 (R 6 ) 11^/11 (5P-3) 
L(3p-1) 

On combine alors toutes les inégalités précédentes pour obtenir (toujours lorsque 

P > Pc) 

7p-9 



T-Lpot{I) — ^ll/ll2i P (R6)ll/llip P (]R6)^cm(/) 2 

ce qui termine la démonstration de la seconde inégalité. □ 

Remarquons que l'énergie potentielle peut être ainsi contrôlée par interpolation entre 
les invariants du système et l'énergie cinétique, avec une puissance strictement inférieure 
à 1 de l'énergie cinétique (plus précisément 1/2 dans notre cas). Par analogie avec les 
équations aux dérivées partielles dispersives et en particulier l'équation de Schrodinger, 
dans le cas d'une telle équation où le hamiltonien est composé de deux termes de signe 
opposé qui peuvent diverger tout en s'équilibrant, on parle d'équation focalisante ; et 
lorsque l'on peut établir un tel contrôle de l'énergie potentielle, on parle d'équation 
sous- critique. Lorsqu'un tel contrôle est possible mais avec une puissance 1 de l'énergie 
potentielle, on parle d'équation critique. Lorsqu'un tel contrôle n'est possible qu'avec 
une puissance strictement plus grande que 1, on parle d'équations sur- critiques, pour 
lesquelles on s'attend en général à des phénomènes d'explosion en temps fini. 

Dans le cas du système de Vlasov- Poisson en dimension 3 qui est sous-critique, 
comme le suggère le résultat que l'on vient d'établir, la difficulté principale pour mon- 
trer des propriétés de stabilité est le contrôle de l'énergie cinétique. Nous renvoyons à 
jLMRll ILMR21 ILMR41 ILMR51 ILMR7] pour une étude systématique de la formation de 
singularités dans des cas critiques (dimension 4 ou dimension 3 relativiste), ainsi que 
pour des commentaires plus précis sur le parallèle avec les équations de Schrodinger. 
Ces modèles critiques n'ont toutefois pas d'interprétation claire au niveau physique. 
Nous nous concentrerons dans cet exposé sur les résultats qui concernent le système 
AH}-© en dimension 3. 

2.4. La théorie de Cauchy 

La théorie de Cauchy est relativement développée pour cette équation non-linéaire. 
C'est même l'un des rares systèmes non-linéaires fondamentaux de la physique pour 
lequel on sait construire des solutions globales et uniques sans hypothèse perturbative 
en dimension 3. Il est facile de se convaincre des points suivants : 

- on souhaite montrer une régularité suffisante sur le champ de force moyen pour pou- 
voir définir des courbes caractéristiques pour les équations différentielles ordinaires 
correspondantes et établir des contrôles sur ces dernières ; 



1044-10 



- cette régularité sur le champ dépend du contrôle sur des moments en vitesse 
f f\v\ k dxdv de la solution, donc en définitive du comportement de la distribu- 
tion pour les grandes vitesses. 

Le contrôle de ces < particules à grandes vitesses » a finalement été obtenu de 
deux manières différentes il y a une vingtaine d'années. D'une part Pfaffelmoser [Pf] a 
découvert comment propager en temps une borne sur la taille du support de la solu- 
tion, ce qui implique le contrôle recherché et permet de construire des solutions globales 
C 1 à support compact. D'autre part et indépendamment, Lions et Perthame |LP] ont 
découvert comment propager en temps des bornes directement sur les moments en vi- 
tesse de toute solution dans L°° et de hamiltonien borné, ce qui leur permet de construire 
des solutions uniques et globales dans L 1 ((l + \v \ k ) àxàv) D L°° pour k suffisamment 
élevé, modulo une hypothèse assez faible de régularité initiale sur la densité spatiale 
p in . Nous renvoyons ensuite aux articles [Ml M IHc^îîl IWo-iil iRôl PJPl iLôl IMMPl ITâ] 
pour les développements ultérieurs de ces théories. 

Par ailleurs, sous des hypothèses plus faibles sur les données initiales 

fin > o, fin e L 1 n l°°(m 6 ), M 2 f in e L 1 ^ 6 ), 

des solutions faibles globales ont été construites par Arsen'ev |Ar] (voir également 
|IN[ IHH] ) ; cependant leur unicité reste ouverte. Ces solutions sont également des solu- 
tions renormalisées au sens de DiPerna-Lions |DLlj. IDL2] et on peut leur associer des 
caractéristiques généralisées au sens de DiPerna-Lions |DL3j . Elles sont continues en 
temps à valeur dans L 1 (M 6 ) et vérifient l'inégalité suivante 

Vt>0, H(f t )<H(f in ) 

(la conservation exacte du hamiltonien est ouverte). Les théorèmes de Lemou, Méhats 
et Raphaël que nous présenterons dans la dernière section ont pour objet ces solutions 
faibles. 

2.5. L'étude des solutions stationnaires 

Nous cherchons maintenant les solutions stationnaires du système OH])-©- H en existe 
en fait une infinité, dont nous ne connaissons et ne savons caractériser de manière 
satisfaisante qu'une petite partie. Une telle solution f° doit vérifier 

Vi,«el 3 , v- VJ° - V x U(x) ■ V v f° = 0, 

Vi G l 3 , A x U = p f o. 

On voit que la première équation inciterait à considérer une fonction /° ne dépendant 
que des invariants du mouvement du système différentiel x"(t) = — V X U, au moins si 
ce dernier est complètement intégrable. Cependant, on ne connaît rien a priori sur ce 
système différentiel, puisqu'il dépend de la solution f° recherchée elle-même. 
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2.5.1. Le théorème de Jeans. - - On sait néanmoins résoudre ce problème dans le cas 
de solutions stationnaires dites à « symétrie sphérique ». Cela correspond à l'invariance 
de la solution stationnaire par les rotations du référentiel, soit 

V (x, v) G M 6 , f°(Ax, Av) = f°(x, v) 

pour toute matrice 3x3 orthogonale A. Cela implique donc que f° = f°(\x\, \v\,x ■ v). 
Les équations du mouvement 

x'(t) = v, 

v'(t) = -V x U = -^U'(\x\) 

possèdent alors toujours (c'est-à-dire indépendamment de /°) un invariant supplé- 
mentaire : 

d , , , , , N , v . , U'(\x\) 

— (x(t) A v(t)) = v(t) A v(t) ^Al x(t) A x(t) = 0. 

Le système est alors complètement intégrable (l'espace des phases possède trois degrés 
de liberté), et l'on peut montrer que les équilibres s'écrivent sous la forme 

( f°(x,v) = f(\x\,\v\,x-v) = F(E,L), 



(12) 



12 



\V\ 

E(x,v) := Efo(x,v) = + 4>fo(t,x) 
L(x, v) := \x A v\ 2 . 



C'est le théorème de Jeans |Jel[ IJe2j appliqué aux systèmes sphériques (voir la dis- 
cussion |BT[ section 4.4]). Un cadre et une démonstration mathématiques rigoureux 
peuvent être consultés dans [BFH] , fondés sur la résolution des équations différentielles 
non-linéaires radiales correspondantes. 

Lorsque F = F(E) ne dépend que de E, on parle de modèles sphériques isotropes, 
et lorsque F = F(E,L), on parle de modèles sphériques anisotropes. Ces termes s'ex- 
pliquent par le fait que la matrice (tenseur) de dispersion en vitesse 

— / ViVjF(E, L) dv, 1 < i, j, < 3 

est un multiple de la matrice identité lorsque F = F(E) alors que, par exemple, ses com- 
posantes diagonales sont différentes lorsque F = F(E, L) dépend de L (cf. à nouveau, 
ici et pour la suite de cette section, \BT\ section 4.4]). 

Dans le cadre de tels systèmes à symétrie sphérique, l'équation fondamentale qui 
gouverne les équilibres de systèmes stellaires est 



(13) 



A x (f) = p -4 ( r W = f F (— + <p, \x A vl 2 } dxdv. 

Jr6 V 2 J 



C'est une équation elliptique que l'on peut voir comme une équation non-linéaire sur le 
potentiel <fi étant donné le profil de l'équilibre sphérique F, ou bien comme une équation 
non-linéaire sur F étant donné 6. 
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2.5.2. Les modèles poly tropiques. - Supposons tout d'abord F = F(E) isotrope, et 
partons de la formule la plus simple pour obtenir un profil F à support compact, c'est- 
à-dire une fonction puissance tronquée : 

(14) /(s, v) = F(E) = C F (E - E) n + , neR, 

où z + := max{z, 0} désigne la partie positive d'un réel z G R, et Cf est une constante 
donnée. Un calcul élémentaire montre que, par résolution de l'équation (fl~3]) . la densité 
spatiale associée est 

p = c n {E - <pyi +3/2 

avec une constante c n finie si n > —1 et la formule simple lorsque n > : 

2 3 / 2 7rn! 
Cn '- CF (n + 3/2)!' 

On voit ainsi que ces modèles ne peuvent jamais être répartis de manière homogène car 
une densité spatiale constante impliquerait n = —3/2. L'équation obtenue alors sur le 
potentiel est l'équation dite de Lane-Emden 

r z ar \ dr / 

On parle de modèles polytropiques car la distribution spatiale est alors la même que 
pour un gaz auto-gravitant vérifiant la loi de pression des gaz polytropes p = este p 7 
avec 7 = 1 + l/(n + 3/2). Le cas n > 7/2 conduit à des solutions de masse infinie et à 
support non compact en ajustant le paramètre E . L'équation sur le potentiel n'admet 
pas de solution simple en terme de fonctions élémentaires en dehors du cas n = —1/2 
(équation linéaire de Helmhotz). Le cas limite n = 7/2 est intéressant car on obtient 
une masse finie mais un support non compact ; c'est le modèle dit de Plummer. 

2.5.3. Le modèle de King. - - D'après les observations, l'énergie est très concentrée au 
centre des galaxies. Cela incite à considérer de grandes valeurs de n dans le modèle 
polytropique. Considérons ainsi tout d'abord formellement le cas limite « n = oo ». 
Bien sûr les équations que nous avons utilisées pour trouver la solution stationnaire 
ci-dessus ne font plus sens, mais si on suit le parallèle avec les gaz auto-gravitants, 
cela doit correspondre au cas d'un gaz vérifiant la loi de pression p = este p 7 avec 
7 = 1, c'est-à-dire un gaz isotherme. En suivant les considérations usuelles de mécanique 
statistique, on obtient des distributions de loi exponentielle. C'est le modèle dit de la 
sphère isotherme. Il est cependant de masse infinie et de support non compact. La 
littérature physique a donc introduit |Kij la troncature suivante 

(15) / = F(E) = (e E °~ E - 1) + 

que l'on appelle modèle de King. Ce modèle reproduit les fortes valeurs d'énergie ob- 
servées pour les petits rayons, tout en étant compatible avec les contraintes de masse et 
support. Ce modèle sera l'un des protagonistes principaux des travaux mathématiques 
récents, car c'est peut-être le modèle sphérique isotrope le plus important d'un point de 
vue physique, et l'étude de sa stabilité échappe aux techniques variationnelles fondées 
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sur les fonctionnelles d'énergie-Casimir. Une des réussites marquantes du programme 
de recherche de Lemou, Méhats et Raphaël est d'obtenir dans [LMR9j la stabilité de 
ces modèles pour des perturbations générales. 

Remarquons ici que l'on voit apparaître en filigrane un phénomème mystérieux, à 
savoir le fait que les équilibres sont analysés par des considérations de mécanique sta- 
tistique collisionnelle, soit irréversible. Comprendre pourquoi cette démarche produit 
des résultats corrects, au moins partiellement, rejoint la question fondamentale de com- 
prendre comment les galaxies se « thermalisent » sans collision. C'est l'objet des théories 
d'amortissement et de relaxation violente initiées par Lynden-Bell |L-Blj. IL-B2] . 

2.5.4- Autres modèles. - - Si l'on se prescrit au départ la densité spatiale et que l'on 
cherche à calculer /, on obtient la formule dite de Eddington |Edj . dont les deux modèles 
célèbres qui en sont issus sont le modèle dit R 1 ^ 4 et V isochrone d'Hénon. Le cas aniso- 
trope est moins clair, mais l'on peut citer le modèle dit d'Osipkov-Merritt qui généralise 
le modèle R 1 / 4 , ainsi que le modèle dit de Michie qui généralise le modèle de King (nous 
renvoyons à |BT| Chapitre 4]). 

2.6. La conjecture de stabilité pour les modèles sphériques 

Nous devons d'abord discuter de la notion de stabilité utilisée. Au vu des invariants 
du système (fonctionnelles de Casimir et hamiltonien macroscopique), il est naturel de 
mesurer la stabilité de la manière suivante : pour tout e > 0, il existe rj > tel que 

(16) V(f m J°) := \\f m - /°|| L1(R6) + \%Uin) ~ H(f°)\ < V 

Vt>0, V(f t ,f°)<e. 

(Les fonctionnelles de Casimir sont correctement capturées par les espaces de Lebesgue 
en terme de distance, et le choix de L 1 (M 6 ) semble arbitraire, mais sera de toute façon 
complété plus loin par une hypothèse de borne L°°(IR 6 ).) 

Cependant, il faut tenir compte du groupe de symétries du système d'évolution et 
voir si ces dernières sont « capturées » par la distance utilisée pour mesurer la stabilité. 
Il s'avère que le système de Vlasov- Poisson flH])- ([9]) admet un large groupe de symétries : 
si f(t, x, v) est solution, alors pour tous t G M, x G M 3 , Ào G fio G 

u v A*o , ( t + t x + x 
g(t, x, v) := -g / -t , — : , Ho v 



Aq V Ao /io Ao 

est encore solution, et on a également l'invariance par translation galiléenne (ce qui ne 
fait que traduire sur ce modèle statistique l'invariance galiléenne correspondante de la 
mécanique classique) : si f(t,x,v) est solution, alors pour tout v$ G M 3 

g(t, x, v) := f{t, x + tv , v + v ) 



est encore solution. Il est clair que la <« distance » utilisée V est également invariante 
par ces transformations lorsque Xq — fj,Q — 1. Nous devons donc reformuler la stabilité 
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de la manière suivante : il existe une fonction de translation z = z(t) G M. 3 telle que 

(17) V(f m J°)< v V*>0, V(f t ,f(--z(t), •))<£. 

C'est la notion de stabilité orbitale. Ajoutons que, dans le cas de perturbations à symétrie 
sphérique, le paramètre z(t) est nécessairement toujours nul. 

Pour les solutions stationnaires à symétrie sphérique, il est très largement discuté 
dans la littérature physique [BTj que la stabilité est liée à la monotonie par rapport à 
l'énergie microscopique E = Efo(x,v). Plus précisément, on peut formuler les conjec- 
tures suivantes : 

Conjecture 2.2. — Les galaxies f° = F(E,L) à symétrie sphérique anisotropes et 
décroissantes (8eF < sur le support de F) sont orbitalement stables par perturbation 
à symétrie sphérique pour le système d'évolution (jHJ)-©- 

Conjecture 2.3. — Les galaxies f° = F(E) à symétrie sphérique isotropes et 
décroissantes (ÔeF = F' < sur le support de F) sont orbitalement stables par 
perturbation générale pour le système d'évolution (JH])-©- 

Les conjectures correspondantes pour le problème linéarisé ont été démontrées dans 
les années 1960 et 1970 à la suite des travaux fondateurs d'Antonov |Anll IAn2j . mais 
leur résolution pour la dynamique non-linéaire correcte vient d'être obtenue par Lemou, 
Méhats et Raphaël respectivement dans les articles |LMR8] et |LMR9j (nous renvoyons 
également à [GL] pour une réponse partielle importante sur la conjecture 12.2p . Nous 
nous concentrerons sur la résolution de la conjecture 12.31 (article |LMR9] ). tandis que 
nous renvoyons à |LMR8] pour la résolution de la conjecture I2.2[ qui implique les mêmes 
idées essentielles. Et lorsque nécessaire, nous illustrerons les méthodes et résultats avec 
les deux modèles principaux : modèles polytropiques f|T4|) et modèle de King ffTol) . 
Ajoutons enfin qu'en l'absence de symétrie sphérique de la solution stationnaire, la 
question de la stabilité est ouverte, même pour le problème linéarisé. 



3. LE PROBLÈME LINÉARISÉ 

3.1. La dynamique linéarisée 

Il est possible d'adopter une vision géométrique de l'espace des solutions du système 
([H])-© en terme de structure de Poisson, avec les invariants définissant une feuille 
symplectique. Nous n'en aurons pas besoin dans cet exposé, et nous renvoyons aux 
références [W£H IYMCI IMPÏI IMP2] et |A"Ml p. 302] pour des approfondissements sur 
cet aspect. Nous mentionnerons uniquement la notion de perturbation dynamiquement 
accessible qui a joué un rôle dans le développement de la théorie perturbative. 
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Si l'on considère une solution stationnaire f°(x, v) on peut linéariser le système (JE])-© 
autour de cette solution / = f° + s h pour obtenir l'équation linéaire suivante 

' d t h + v- W x h - V x (f) f o ■ V„/i = W x <p h ■ VtJ , t E R+, x G M 3 , v E M 3 , 



Il est aisé de construire des solutions régulières uniques et résoudre le problème de 
Cauchy pour ce système d'équations aux dérivées partielles linéaires, par des méthodes 
de point fixe et en utilisant l'inégalité d'interpolation fondamentale de la proposition ^. Il 

On peut réécrire ce système d'évolution de manière plus compacte en 



Il est important de remarquer que la fluctuation h n'est en général pas positive 
et que la structure hamiltonienne de l'équation (jHJ) n'a pas survécu au processus de 
linéarisation. Le système linéarisé (I18p en hérite toutefois la propriété suivante. Si l'on 
définit une perturbation initiale dynamiquement accessible comme étant créée à partir 
de f° par un hamiltonien extérieur donné H ext , on obtient formellement 

/ = f + e h + 0(e 2 ) et H = E f o+sH ext avec h\ t=0 = 0. 
Le premier ordre de linéarisation en e de l'équation sur h donne 



rôle de la direction de dérivation dans le processus de linéarisation. On montre alors 
dans la proposition suivante que la forme h t = {gt, f } est préservée au cours du temps 
pour l'équation linéarisée ffT8]) sur h, soit que l'on peut restreindre l'équation linéarisée 
à l'« espace tangent » des conservations hamiltoniennes, i.e., de la feuille symplectique. 

Proposition 3.1. - - Si g est une solution (régulière) de /'équation génératrice sui- 
vante 



(18) 




d t h + {h,E f o} = V x (p h ■ 
h(t = 0,x,v) = ho(x,v). 




(19) 





alors h = {g, f } est une solution régulière de f[TÏÏ|) . 



Remarquons dès lors, ce qui nous sera utile par la suite, qu'il est équivalent de chercher 
une solution sous la forme h = {g, f } ou sous la forme h = {g, Ep} : il suffit en effet 
de multiplier g dans la première forme par l/F'(Ep) pour obtenir la deuxième. 
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Démonstration. — On calcule l'équation d'évolution sur ht : = {gt, f } '■ 

dth = {d t g, f} = -{{g, E f o}, f} + {<f> h , f} 

= {{E f o,f }, g} + {{f, g}, E f o} + f} 

où l'on a utilisé l'identité suivante sur le crochet de Poisson 

{{A, B}, C} + {{B, C}, A} + {{C, A, }, B} = 0. 

On utilise alors {Efo, f } = (par définition d'une solution stationnaire), et {0^, /°} = 
V x 0ft • V v f°, ce qui permet de conclure à 

dth + {h 7 Efo} = V x <f> h -V v f 
qui est bien l'équation recherchée, avec la donnée initiale hi n = {gi n , f }- D 

3.2. Energie libre du problème linéarisé 

Dans le cas de modèles sphériques, l'article |KO] a pour la première fois proposé 
une énergie libre préservée au cours de l'évolution par le système de Vlasov-Poisson 
linéarisé dans le cas d'interactions de Coulomb, pour les plasmas. Cette fonctionnelle 
sera ensuite adaptée au cas gravitationnel du système ([18]) par Antonov |Anl] IAn2j , et 
ce dernier l'utilisera pour formuler un critère de stabilité. Nous ne présentons ici que 
le cas isotrope f° = F(E) pour simplifier l'exposition et nous rappelons la notation 
E = Efo(x, v). 

Proposition 3.2. — La quantité 

F{h) = -f JMLdzdu- / |V^| 2 dx 



est conservée par les solutions régulières du système (jl8j) dont le support est inclus dans 
le support de F' . 

Le quotient \h\ 2 /F'(E) doit ici être compris comme étant égal à zéro lorsque les 
numérateurs et dénominateurs s'annulent. 

Démonstration. — On souhaite différencier en temps l'expression de F {h). On calcule 
dx dv = — 2 / - ' ; — dx dv + 2 f — '^t^.^J — dx dv 



dt 7k 3 xr 3 F'(E) J R 6 F'(E) J R 6 F'(E) 



puis on remarque que 



-2 / dxdv + 2 I v-V x (j) h hdxdv 



-^ 2 {E,E} = 



F'iE)' J F'(E) 
et ainsi par intégration par parties 

o f i h ' E } h dxdv= f ^;^} dxdv = - I h 2 {——^E}dxdv = 



F'(E) UF'(E) J R6 \F'(E) 
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d'où 

A f _M_ 

di 7r3 xR 3 F'(E) 
On calcule par ailleurs 



dxdv = 2 / -y • V x .0h /idxdu. 



f d f d 

\V x <f> h \ 2 dx = 2 V x .— <j) h ■ V x (p h dx = -2 ^-A^ </» h dx 

CirjRS JR3 Ot jR3Ct 

= -2 [ ^-ph4>hdx = -2 [ ( [ T^hdv J ç^dxdf 

— —2 I v ■ V x hcf) h dx = 2 v ■ V x 4>hhdx. 
La somme de ces deux dernières équations fournit le résultat. □ 

3.3. Interprétation variationnelle de l'énergie libre 

Quelle est l'origine de cette énergie libre et surtout quel est son lien avec le hamiltonien 
du problème non-linéaire ? Si l'on compare le hamiltonien en deux fonctions différentes 
f° et / on obtient par un calcul élémentaire la formule 

(20) H(f) = H(f°)+ [ fÇ + ^oV/-/°)dxdt;-i/ |V^,-V,^o| 2 dx. 

JR6 V 1 J l J9? 

En effet 

+ J r6 (^Y + (f - /°) dxdv + fy> f f - <j> p {f - /°) - ^ /0 /°) dxdv 
= H(f°)+ [ E(f-f°) dxdv- f (hv^ff-V^fo-V^f + hv^foA dxdv 

JR6 Jr3 z J 

d'où le résultat. 

Il est clair cependant qu'aucune fonction f° ne peut constituer un point critique de % 
à cause de la partie linéaire en (/ — f°) dans la formule (|2"0"j) ci-dessus, qui n'est jamais 
l'application nulle. Cela s'explique par les invariants du système. Pour y remédier, 
une idée introduite par Arnold |Arl} IAr2} IAr3j dans les années 1960 dans le cas de 
l'équation d'Euler incompressible bidimensionnelle, est de considérer une fonctionnelle 
dite d'énergie- Casimir qui combine le hamiltonien et une fonctionnelle de Casimir bien 
choisie. Cette idée resurgira plus loin lorsque nous aborderons la stabilité non-linéaire. 

On se donne alors une fonctionnelle d'énergie-Casimir générale définie par une fonc- 
tion C : R + ->• R + avec C(0) = : 

Uc{f):=U{f)+ [ C(f) dxdv, 
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et l'on cherche à ajuster la fonction C afin que cette fonctionnelle admette bien f° 
comme point critique. Par un calcul élémentaire et en développant Tic autour de f° on 
obtient 

n c (f)=Uc(f )+ [ O-^ + tfo+C'if^if-f^dxdv-l f \V x <j> f -V x ^ P \ 2 àx 

+ U C"{f){f-ffdxdv + ... 

=: Uc{f) + DU c {f)[f - f] + \D 2 Uc{f)[f - f] + • • • 



\vy 

T 



On voit alors que le choix formel 

C\f) = —E = - 
permet d'annuler la partie linéaire et d'obtenir exactement 

D 2 Hc(f )[f-f } = Hf-f )- 
En effet puisque f° = F(E) ne dépend que de E on a 

C'(F(E)) = —E =>- C"(F(E))F'(E) = -1 soit C"{F{E)) = -j^- y 

On peut donc interpréter cette énergie libre comme la hessienne de la fonctionnelle 
d'énergie-Casimir au point f°, pour une contrainte C correctement choisie en fonction 
de /°. 

3.4. Première approche naïve de la stabilité linéaire par interpolation 

Donnons tout d'abord une approche « naïve » de la stabilité qui rappelle l'inégalité 
d'interpolation de la proposition 12.11 pour l'équation non-linéaire. On peut en effet 
interpréter à nouveau la fonctionnelle T comme une énergie que l'on décompose en 
deux termes positifs : 

J={h) = F cin (h) - Fpotih) avec 
Fcin(h) := - I — — dxdv, FpotW := / \V x (j) h \ 2 dx 

JR3 X R3 P {&) JR3 

et l'on cherche à contrôler la partie J-p ût à partir de Tcin dans une estimation sous- 
critique. La proposition suivante est inspirée de |BMR] . 

Proposition 3.3 ( |BMRj ). - - Si la solution stationnaire f° vérifie 

(21) / F'(E) dxdv < +oo 

alors 

et l'on peut démontrer la stabilité dans un espace de solutions telles que ph G L°°(IR 3 ). 
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La preuve de l'inégalité fonctionnelle est immédiate, et le résultat de stabilité linéaire 
en découle une fois les solutions correctement construites. Malheureusement l'hypothèse 
(12T]) ne permet pas de traiter les modèles physiques les plus intéressants. Il va donc falloir 
se tourner vers une estimation globale de T sans décomposition; autrement dit, une 
estimation de convexité sur la fonctionnelle d'énergie-Casimir T-Lq- 

3.5. L'inégalité de coercitivité d'Antonov 

Antonov |AnltrA"ïï2] propose alors de considérer la fonctionnelle J 7 dans son ensemble, 
et de montrer qu 'elle est positive, une fois restreinte aux perturbations dynamiquement 
accessibles h = {g,f }. Pour être exact, l'argument (formel) original d'Antonov n'est 
pas fondé sur les perturbations dynamiquement accessibles, mais sur une décomposition 
de la perturbation en partie paire et partie impaire par rapport à la variable de vitesse. 
Nous renvoyons également aux travaux ultérieurs dans la littérature physique qui ont 
développé l'idée d'Antonov jDFBl iDFl ISdFLRPl [KSI IFA] . 

L'argument historique d'Antonov consiste à considérer une solution h du problème 
linéarisé (fT9~|) puis à la décomposer en partie paire et partie impaire par rapport à la 
variable de vitesse v , soit h = hi + h 2 avec 

hi(x, v) = - (h(t, x, v) + h(t, x, -v)) , h 2 (x, v) = - (h(t, x, v) - h(t, x, -v)) . 

On obtient ainsi formellement les deux équations couplées suivantes 

dth% + {h 2 , E} = 0, 

d t h 2 + {h 1 ,E} = V x (j )hl -V v f . 

Puis en remplaçant le terme h\ dans la deuxième équation, on déduit V équation de 
pulsation suivante sur h 2 

(22) d 2 tt h 2 = {{h 2 , E}, E} - V x (j){h 2 , E } ■ V,/ . 

Par analogie avec l'équation différentielle y" = Xy, on voit que, si l'opérateur du membre 
de droite est négatif, on espère obtenir des solutions oscillantes bornées (ce qui justifie 
le terme de « pulsation »). 

On calcule alors l'énergie de ce système. 

Proposition 3.4. - - On suppose E' < sur son domaine, et h g une solution 
de (122|) à support inclus dans celui de F' . Alors h vérifie 

soit 

Tt{-LjW) Axiv + H{h ' f)) )^ 
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Démonstration. — En intégrant contre dth/F'(E) l'équation de pulsation on obtient 

/ d%h — \ \-. dx dv — -, — / [_,.*/ In dx dv 
J R e tt F'(E) dt2j R6 F'(E) 

pour le premier terme, puis 

J., 1 F'(E) 2di7 R , 

pour le deuxième terme où l'on a utilisé une intégration par parties et le fait que 
{1/F'(E),E} = 0, et enfin pour le dernier terme 

- I Vx<j>{h,m ■ V p /° * dxdv = - / {0 {/ïie} ,£ , }<9i/idxdî; 

= / ^{h,s}{^/i,J5}da;dt; = / <f){ h ,E}d t p{h,E} dx 

= - [ Vx<fi{h,E} ■ Vxdt<fi{h,E} dx = ~-tt [ I ^x<f>{h,E}\ 2 dx. 
On obtient le résultat souhaité en combinant les trois précédentes égalités. □ 

On voit donc que, dès lors que l'on sait montrer la positivité de l'énergie libre sur les 
perturbations admissibles {h 2) f } issues d'une fonction /i 2 impaire, on peut en déduire 
l'impossibilité d'instabilités sur la composante impaire h 2 puisque l'on a alors la somme 
de deux termes positifs qui est conservée au cours du temps. C'est le critère de stabilité 
d'Antonov |Anlj . Le reste de l'argument de |Anlj est moins clair en ce qui concerne 
le contrôle de la composante paire h\, et semble conditionné à certaines hypothèses ad 
hoc pour éviter la possibilité d'instabilités avec une croissance polynomiale en temps. 
Dans l'article suivant |An2j . Antonov donne une preuve de cette propriété cruciale de 
positivité de J-"({/i2, E}) avec h 2 impaire, en résolvant le problème de minimisation 
associé. 

Nous allons maintenant présenter les analyses mathématiques récentes de ces ques- 
tions. Tout d'abord, nous présentons un énoncé et une preuve élégante issue de |GR3] de 
la coercitivité découverte par Antonov, i.e., la positivité de l'énergie libre sur les fonc- 
tions de la forme {h 2 , f } avec h 2 impaire. Du fait de l'importance de cette propriété 
d'un point de vue historique, ainsi que pour la genèse des résultats dont il est question 
ici, nous présentons une preuve détaillée. Une autre approche de cette propriété et de 
sa preuve, développée dans |LMR8t ILMR9] . sera discutée dans la section m 

Proposition 3.5 ( |GR3j ). - - On suppose F' < sur le domaine de f°, et h e à 
support inclus dans celui de F' , à symétrie sphérique et impaire par rapport à la variable 
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de vitesse. Alors on a 

r ( „ r h ï 2 «r), 

1 f h 2 dxdv > 



(23) H{Kf°})>-[ F'{E)[ 



x ■ v\ 



X ■ V 

où r = \x\ et 4>'p(r) désigne la dérivée radiale de = cj)p{r). 

Démonstration. — Observons tout d'abord que, pour un modèle sphérique f° = F(E), 
les fonctions <pp et pp sont radialement symétriques, et on a 

A x .0/o = pp =^ — [r 2 <p'p(r)]' = pp(r) 
d'où, en intégrant de à r, on obtient 

1 f r 

Vr > 0, (j>' f o(r) = — rîp f o(r*) dr* > 
r Jo 

et le potentiel créé par la solution stationnaire f° est croissant. On déduit donc 
immédiatement la positivité du membre de droite dans la proposition. 
Rappelons la formule de l'énergie libre 

F({hJ°}):=- [ l{k Jl\ 12 dxdv - [ \V x <f> {h>n \ 2 dx. 

Nous allons tout d'abord contrôler par au-dessus l'opposé du second terme. On observe 
d'une part que la symétrie sphérique de h implique immédiatement celle de cj){hj }, et 
d'autre part que 

/ {h,f°}dv= [ (V x h ■ V v f° - V v h ■ VJ°) dv 

= Va ■ ( [ hV v f°dv) = V x ■ ( f vF'(E)hdv) . 

v^r 3 y v^r 3 y 

Par conséquence, en utilisant à nouveau A x (p = r _2 [r 2 0']', on obtient 

4>'{ h ,m( r ) = ~2 I (S {h(x*,v),f(x*,v)}dv ) r 2 dr* 
r Jo Vr 3 y 

= ^L v ' ' (b F ' (E)hAv ) àx = ^LL fëf) W>*>***) 

où du désigne la mesure uniforme sur la sphère, et dans la première ligne = r*. 
Puis, l'intégrande est à nouveau invariant par rotation sur la variable x et on en déduit 

V W = r, <j>' {hJQ} {r) = [ (^f) F'(E)hdv. 

Jm.3 v fi y 

On aboutit au contrôle suivant sur le second terme de l'énergie libre par inégalité de 
Cauchy-Schwarz : 

f |V^ {hJ o } | 2 dx< / (-/ (^-Y F'(E)dv) (-[ F'(E)\h\ 2 dv) dx. 

jr 3 7r 3 v >^r 3 v \ x \ y y v Vr 3 y 
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Étudions plus précisément la première intégrale en v du membre de droite. En chaque 
point x, on décompose orthogonalement la variable v en 



x ■ v \ x ( x A v 



x 





M + 






\x\ 


/ \x\ 




/ " M 



et 




-2 / / ^F' 



«H , Kl 2 



2 2 

zWr) I dvu dvi 



o dvu dv i 




2 2 



I ?r + ^4- + 0/o (r) ] dvu dv± = I f°(x,v) dxdv = p P {r). 



2 2 



On a donc 



/ \V x (t> [hJ o } \ 2 dx<-( F'(E)p f o(r)\h\ 2 dxdv. 



On calcule maintenant le carré du crochet de Poisson suivant 



|{*,Al 2 

= (x • v ) 2 



(x ■ v) 



(x ■ v) 
(x ■ v) 2 



2 / , \ 2 
Ai 



[x ■ v) 

2 



\{(x-v),f}\' + 2h 



[x ■ v) 

h 



(x ■ v) 



[x ■ v) 



(x ■ v) 
(x-v){(x-v),f},f° 

h 2 



J°){(x-v),f} 



[x ■ v) 



{{(*■«), A, /°} 
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ce qui donne, lorsque l'on intègre contre 1/F'(E) : 



\{h,f°}\ 2 
F'(E) 



dx dv 



+ 



(x ■ v) 2 
F'(E) 



h 



[x ■ v) 



f 



dx dv 



F'(E) I \(x-v 



h 



Ir 



,F'(E) \(x-v 
f (x-v) 2 F>(E) l 



(x-v){(x-v)J°}J 0[ j dxdv 

{{(x-v)J°}J°}dxdv 
h 



E 



dx dv 



+ 



h 



[x ■ v) 



(x ■ v) ' 

(x ■ v ){(x ■ v), f }, E } dx df 
h 2 



[ F\E) (-^--]{{(x-v),E},E}dxdv. 
Jm. 6 \ \ x ' v ) J 



Le deuxième terme du membre de droite s'annule par intégration par parties, et comme 
{{(x ■ v), E}, E} = -(x ■ v)^ (r) - 3^% (r) 

= Wo(r)r 2 + 20' /o (r)r) - ^% (r) 



(r 2 c/)/o(r)) <j>' f0 (r) = -(x ■ v)p f o(r) 0/o(r), 



on obtient pour le troisième terme 
h 2 



F'{E) (j^—A {{(x-v),E},E}dxdv= [ F'{E)h 2 (p f0 + < ^-\ dxdv 
V \ x ' v ) J Jr6 \ r J 



et donc 



dx dv 



On conclut la preuve en combinant cette égalité avec l'inégalité précédente sur 

J u3 \V x (j) {h ,fO } \ 2 dx. □ 



3.6. Une inégalité de coercitivité précisée à la Weinstein pour les polytropes 

Pour conclure cette section, nous allons maintenant considérer le modèle polytropique 
(fT4"|) et présenter l'analyse linéarisée rigoureuse effectuée dans |LMR3j . Nous renvoyons 
à la section suivante pour la caractérisation variationnelle des polytropes de |LMR4j qui 
est utilisée pour montrer la positivité au sens large de la fonctionnelle d'énergie libre 
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du problème linéarisé, et nous montrons comment les auteurs de |LMR3] en déduisent 
une inégalité de coercitivité précisée. 

Nous présentons ce résultat et une ébauche de preuve car cette dernière contient 
l'une des idées de la méthode non- linéaire de |LMR9] : 1' inégalité de coercitivité sur 
l'opérateur de Schrodinger (128]) ci-dessous, qui sera utilisée pour pouvoir traiter les 
perturbations non radiales dans le cas non-linéaire. Ce travail est inspiré de l'étude 
par Weinstein |We-ii-14 IWe-ii-2] dans les années 1980 de la stabilité des solitons pour 
l'équation de Schrodinger. 

On considère une solution stationnaire polytropique ( IT4"|) avec n = — 1) et 
E = -1 : 

/ i..|2 \V(P-1) 

(24) f = F(E) = (-1 - Ef/^ = (-1 - a- - ^oj 

pour p > p c = 9/7 supérieur strictement à l'exposant critique que nous avons déjà ren- 
contré dans l'inégalité d'interpolation de la proposition 12.1} ce qui redonne exactement 
la condition n < 7/2. On définit l'opérateur 

Mh= {-JW) +< ") lK 

restreint au support K de f°. Il est clair que 

{Mh, fc)L2(R3 X R3) = T{h) 

et que cet opérateur linéaire est symétrique dans L 2 (K, d/i) avec la mesure de référence 
fi = 1/\F'(E)\. On peut donc reformuler la question de contrôler par en-dessous l'énergie 
libre en un contrôle de coercitivité sur l'opérateur A4. 

Tout d'abord, nous admettons ici la positivité (au sens large) A4 > de cet opérateur 
pour des perturbations qui ne modifient pas le hamiltonien : 

h<EL 2 (K,dfi) avec f hE dxdv = f h ( — - — h </>/<> ) dxdv = 0. 

Cette positivité découle de la caractérisation variationnelle de la proposition 14.11 par 
saturation d'inégalité de Sobolev de la section suivante. 

Une fois cette positivité acquise, en s'inspirant de |We-ii-lj . Lemou, Méhats et Ra- 
phaël quantifient la coercitivité de l'énergie libre de la manière suivante. 

Proposition 3.6. - - Pour p c < p < +oc, la forme quadratique h y (A4 h, h) est 
continue et auto-adjointe sur L 2 (K, dfi), et il existe une constante ô ne dépendant que 
de p telle que 

(25) (MKk)>_sf K ^-\ 

On retrouve le résultat précédent (avec cependant une constante non constructive) 
lorsque h = {g, E} et g impaire en v, mais l'analyse du défaut de coercitivité pour les 
perturbations qui ne sont pas sous cette forme est ici plus précise. 



hE dx dv 



i=l 



X; h dx dv 
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Ebauche de preuve. — On va raisonner par l'absurde. On définit Z l'ensemble des fonc- 
tions h G L 2 (K, du) telles que 

hE dxdv= / Xihdxdv = / x 2 hdxdv= / x 3 hdxdv = 0. 

JKP JKP Jr 6 

On sait par la caractérisation variationnelle de la proposition U7T] que (Ai h, h) > pour 
h G Z. On suppose alors que 

I : = inf (Mh,h) = 0. 

heZ, \\h\\t2( K ,dri =1 

Étape 1 - Construction d'un minimiseur. On considère une suite minimisante h n G Z. 
Par compacité faible, quitte à extraire, la suite h n converge faiblement vers une fonction 
dans L 2 (K, d/i). Par ellipticité de l'équation de Poisson, on a la convergence forte 
de V x (j) hn dans L 2 (M 3 ). De par la normalisation, on a 

(Mh n ,h n ) = 1 - HV^JI^s) ->■ 

d'où || Vx^h.^ ||l 2 (ir 3 ) = 1 et la limite /loo n'est pas nulle. 

Étape 2 - Le minimiseur est dans le noyau de Ai. Par la méthode des multiplicateurs 
de Lagrange, on déduit du problème de minimisation sous contraintes ci-dessus que 

Aihco G Vect | ^^y x l 1 K,x 2 lK,x 3 l K ,El K ^ . 

En considérant les intégrations de .M/ioo contre successivement hoo, d Xl f°, d X2 f°, d X3 f° 
et x ■ V v f° — 2v ■ V„/°, on annule chacun des coefficients dans la décomposition selon 
cette famille vectorielle, et on déduit finalement que Aihoo = 0. 
Etape 3 - Etude du noyau de Ai. On montre que 

KerM = Vect {d x J°, d x J°, d x .J } . 

L'inclusion de ces trois vecteurs dans le noyau de Ai provient simplement de la 
dérivation par rapport à x%, x 2 et x 3 de l'équation définissant la solution stationnaire 



(26) 



= d Xi ((/°r 1 + t£ + + l K = (j^dxif + 0^/o) 1k- 



Pour l'inclusion réciproque, on considère l'équation 

h 



( 27 ) p l{E) r» 

sur K, et l'on en déduit l'équation réduite suivante sur le potentiel par intégration en 
vitesse 

A x (j)h = Ph= / hdxdv= / F'(E)(j) h dxdv. 
Si l'on définit le potentiel effectif 



V f o :=- [ F'{E) du, 
Jm 3 
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on fait donc apparaître une équation de Schrôdinger stationnaire 
(28) Mh •= (A« + Vfo) <f> h = 0. 

Le noyau de cet opérateur A est alors étudié de manière fine par décomposition selon 
les harmoniques sphériques, avec pour résultat 

Ker(^4) = Vect {d Xl 4>p,d X2 4>p,d X3 (j)p} . 

On déduit donc que (fih G Vect {d Xl (j)fo,d X2 (j)fo,d X3 (j)fo}, puis en utilisant fl27|) et (1261) 
quehe Vect{d x J ,d x J°,d X3 f }. 

Étape 4 - Conclusion. Puisque h G Z, en intégrant la relation linéaire h G 
Vect {d Xl f°, d X2 f°, d X3 f } contre x±, %2 et x$, on obtient successivement que tous 
les coefficients de la combinaison linéaire sont nuls, soit h = 0, ce qui aboutit à une 
contradiction. □ 

Mentionnons qu'au moyen de cette inégalité de coercitivité, Lemou, Méhats et Ra- 
phaël démontrent ensuite dans |LMR3] un théorème de stabilité linéarisée qui énonce 
que, pour toute donnée initiale dans un « espace d'énergie » (correspondant au problème 
de minimisation non-linéaire), le semi-groupe linéarisé croît au plus en 0(t 2 ). Ils donnent 
également une décomposition de l'espace L 2 (K,dfi) qui localise plus précisément les 
modes de croissance algébrique. 



4. LA STABILITE NON-LINEAIRE 

Nous allons maintenant suivre le cheminement des différents travaux sur la stabilité 
non-linéaire. En dehors du dernier travail |LMR9] que nous détaillerons, nous donnons 
seulement les étapes principales. 

4.1. Premières approches variationnelles 

Le grand succès de ce programme conduit par différents groupes indépendants est 
la preuve de la stabilité de tous les modèles polytropiques discutés précédemment. Les 
ingrédients communs à ces différentes approches sont : 

(a) la caractérisation de la solution stationnaire f° étudiée comme un état fondamental 
(< ground state ») d'un problème de minimisation de la forme 

(29) min U{f) 

V(/)=cstes 

pour une certaine fonctionnelle IÂ et un ensemble de contraintes V, qui sont 
préservées par l'évolution non-linéaire; 

(b) la preuve d'une propriété de séparation des états fondamentaux : la solution 
stationnaire f° est isolée parmi les minimiseurs du problème de minimisation 
précédent ; 
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(c) la compacité des suites minimisantes, qui repose souvent sur une technique de 
concentration-compacité [Lil|. ILi2j . sachant que la compacité doit être obtenue 
dans un sens assez fort pour permettre d'obtenir à la limite une solution station- 
naire du système ; 

(d) la preuve de la stabilité est alors fondée sur un raisonnement par l'absurde : on 
considère une suite minimisante qui ne reste pas proche de la solution station- 
naire f°, puis en appliquant (c) on aboutit, à la limite, à un état stationnaire qui 
minimise le problème (|29|) mais qui est différent de f°, ce qui contredit (a)-(b). 

L'idée d'utiliser une fonctionnelle d'énergie-Casimir bien choisie et d'étudier ses points 
critiques et sa convexité a été introduite pour les équations d'Euler incompressibles en 
dimension 2 par Arnold |Arl|. IAr2j. IAr3] , et elle a ensuite été appliquée avec succès aux 
équations de Vlasov-Poisson pour les plasmas dans jHMRWt IRel] (voir également les 
références incluses dans jHMRW] pour les travaux antérieurs de physique sur la stabilité 
formelle pour les plasmas, ainsi que |Gul[ IGu2j pour des travaux mathématiques dans 
le cas de plasmas magnétiques). 

En ce qui concerne le système de Vlasov-Poisson gravitationnel, le premier travail 
précurseur en ce sens est dû à Wolansky |Wo-ij : ce dernier caractérise l'équilibre comme 
le minimiseur d'une fonctionnelle d'énergie-Casimir 



Cette approche a été développée de manière systématique par Guo et Rein |Gu3[ IGR1[ 
IRe2j et a permis d'obtenir la stabilité des polytropes f fT4|) pour < n < 3/2. Les cas 
3/2 < n < 7/2 ont ensuite été étudiés dans jCu4l ICR21 IRCl ISc2l IHâ] . En particulier 
les articles |Gu4[ IGR2] modifient le problème variationnel de la façon suivante : 



Parallèlement et de manière légèrement différente à ces travaux, Dolbeault, Sânchez 
et Soler |DSS] introduisent ensuite un problème de minimisation différent du type 

mmH(f) sous les contraintes / > 0, ||/||i,i(R6) = M, ||/|| L <»( R 6) < 1. 

Ce problème de minimisation se réécrit de manière équivalente et naturelle par satu- 
ration d'inégalité fonctionnelle de type Poincaré reliant l'énergie potentielle, l'énergie 
cinétique, et les normes utilisées : 



Cette approche a permis de traiter le cas formel limite n = dans f lT4|) (en plus de cas de 
polytropes anisotropes que nous n'évoquons pas ici). Sânchez et Soler |SSj ont ensuite 
généralisé cette approche à un espace de contrainte ||/||z,i(R6) = M et ||/||lp(r6) < 1, 
et ont pu montrer la stabilité au sens de la distance L 1 (M 6 ) pour les polytropes ffT4"j) 




min"H(/) sous les contraintes 




mm 




sous les contraintes / > 0, ||/||z,i(R6) = M, \\f\\ L °°(W>) < 1- 
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avec < n < 7/2. Un des apports de ces travaux semble conceptuel : montrer que le 
problème variationnel sous-jacent est relié à des inégalités de type Sobolev optimales. 

Simultanément, Lemou, Méhats et Raphaël [LMRll ILMR4j caractérisent les poly- 
tropes à partir d'une inégalité de type Sobolev optimale correspondant aux inégalités 
d'interpolation de l'équation non-linéaire. Le problème de minimisation en terme de 
fonctionnelle d'énergie et d'espace de contraintes est équivalent à celui considéré par 
Sânchez et Soler. Mais ils font ainsi le lien avec l'inégalité d'interpolation de la proposi- 
tion [2IH et ils effectuent aussi un retour conceptuel à la méthode originelle de Cazenave 
et Lions jCL] pour l'étude de la stabilité des solitons par concentration-compacité pour 
l'équation de Schrodinger. 

Ils démontrent la proposition suivante, que nous avons déjà évoquée et utilisée dans 
l'étude linéarisée pour la preuve de la proposition 13.61 

Proposition 4.1 ( [LMRll ILMR4] ). — Soient p e}p c , +oo[ et f° un polytrope défini 
par (j24|) . Alors le problème de minimisation 

\\\ v \ 2 f\\%(RS)\\f\\%(R6)\\f\\%(Re) „ 1 „ P „ (7p-9) 

mm — — —75 avec t>i = -, u 2 = —, r, 03 = —, r 

î^m IIV.0,111™ 2' 2 3(p-l)' 6(p-l) 



(le dernier coefficient est bien positif du fait que p > p c = 9/7) est atteint sur la famille 
à quatre paramètres 

if ( tJiv^, 7 e m;, A6i;,/i6 m;, x e m 3 . 

En procédant selon les grandes lignes de la stratégie décrite plus haut, les auteurs 
démontrent ensuite dans |LMR4j la stabilité des polytropes (TT41) pour < n < 7/2. 

Cette approche variationnelle a pu traiter de manière satisfaisante les modèles po- 
lytropiques. Elle semblait cependant impuissante à traiter des modèles plus généraux, 
et en particulier le modèle de King. La difficulté est que pour les types de problèmes 
de minimisation que l'on vient d'énumérer, la propriété de séparation des états fonda- 
mentaux (b) décrite plus haut n'est en général plus vérifiée, sans même parler de la 
propriété (c) de compacité des suites minimisantes. 

4.2. Approche directe non-variationnelle par linéarisation 

Il existe essentiellement deux manières d'aborder la question de la stabilité d'un 
système d'évolution non-linéaire : d'une part l'approche variationnelle où l'on exprime 
la solution stationnaire comme solution d'un problème de minimisation qui est invariant 
le long de l'évolution, et d'autre part l'approche « directe » par linéarisation et contrôle 
du reste. Nous allons maintenant parler des travaux s'inscrivant dans cette dernière 
approche. 



4. Nous renvoyons également au travail jCSSj qui étudie selon une stratégie proche les propriétés 
de stabilité orbitale pour l'équation de Nordstrom- Vlasov dans un cadre relativistc. 
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Dans le cas de cette approche directe, il faut tout d'abord quantifier précisément les 
propriétés de stabilité du système linéarisé. L'inégalité de coercitivité d'Antonov ( Î23î> 
est un point de départ naturel pour cela. Cependant il faut ensuite surmonter deux 
difficultés importantes : 

- il faut contrôler les termes d'ordre supérieur ou égal à trois dans le développement 
de Taylor de la fonctionnelle d'énergie-Casimir au voisinage de la solution station- 
naire f° considérée ; 

- et l'autre difficulté est que l'inégalité de coercitivité d'Antonov ( 123]) n'est valide 
que pour les perturbations dynamiquement accessibles de la forme h = {g, f }, et 
que l'on souhaiterait s'affranchir de cette restriction. 

La première tentative d'utiliser cette approche remonte à Wan |Walj mais la preuve 
semble incomplète. La deuxième tentative du même auteur |Wa2j est plus aboutie mais 
semble reposer sur une hypothèse non réaliste de positivité de la fonctionnelle T pour 
toute perturbation h qui exclut la plupart des modèles physiques. 

Le premier article traitant du modèle de King, et suivant cette approche directe, est 
dû à Guo et Rein |GR3j . Les ingrédients clés de ce travail sont : 

- la définition d'une classe de perturbation à symétrie sphérique 

S p = {/ g L X (1R 6 ), / = f(E,L) > tel que 

VC G C 2 (R 2 + ) avec C(0, L) = diC(0, L) = 0, d\C borné, alors 

( C(f,L)dxdv= [ C(f°,L) dxdv) 

qui est : d'une part stable par le système d'évolution non-linéaire du fait que L est 
conservé le long des trajectoires et donc les fonctionnelles f R6 C(f, L) dxdv sont des 
invariants du système et, d'autre part, incluse dans les perturbations de la forme 
h = {g, f } (ce dernier point est démontré en résolvant le système différentiel 
ordinaire associé) ; 

- la démonstration par contradiction d'une propriété de convexité stricte de la fonc- 
tionnelle d'énergie-Casimir "H au voisinage de f° : 

n{f) - n(f°) > CoWv^f - v x <p f 4h m 

pour une constante Co > et avec 

d(f, f°) := n(f) - H(f) + || V,0/ - V x (f> f o\\ 2 LHR3) assez petit. 

Il y a deux limitations importantes dans les résultats ainsi obtenus. D'une part les 
perturbations considérées sont restreintes à la classe Sfo, qui est « trop petite », et en 
particulier incluses dans l'ensemble des fonctions équimesurables à f°. D'autre part, on 
aimerait bien sûr s'affranchir totalement de la contrainte de symétrie sphérique pour 
ces perturbations. 
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La première de ces limitations a ensuite été levée dans le travail |GL] de Guo et Lin. 
Ils démontrent ainsi la stabilité du modèle de King par petite perturbation à symétrie 
sphérique. Les éléments principaux de leur travail sont les suivants. 

- Ils considèrent la fonctionnelle d'énergie-Casimir 

Hc(f) = H(f) + / C(f) dx dv avec C(f) := (1 + /) ln(l + /) - 1 - / 
et la distance associée 

d c (L f°) ■= Hc(f) - n c (f°) + \\v x <f> f - v x <t> f0 ||| 2(M3) . 

- Afin de montrer la coercitivité de cette fonctionnelle d'énergie-Casimir au voisinage 
de f° sans faire apparaître de termes d'ordre supérieur, ils font appel à une inégalité 
de dualité convexe élémentaire mais astucieusement utilisée, qui permet de contrôler 
par en-dessous 

U c {f)-Uc{f) >cste ( [ |V x .0 / _ / o| 2 dx+ / F'(E)\(f) f _ f0 -V(f>f- f o\ 2 dxdv] 

- I \V<p f -p\ 2 dx - 

JR 3 

où les trois points désignent des termes contrôlables par les invariants du système, 
et l'opérateur V est le projecteur sur le noyau de D := v ■ V x — V x 0/o • (c'est un 
opérateur de moyennisation sur chacun des tores invariants du flot complètement 
intégrable associé à cet opérateur de transport). 
-Malheureusement le terme négatif —f R3 \V(fif-fo\ 2 dx dans l'équation ci-dessus 
semble difficilement contrôlable, aussi les auteurs ont-ils l'idée d'approcher 
l'opérateur V par une suite bien construite d'opérateurs Vg de rang fini, et de 
remplacer V par Vt dans l'argument ci-dessus. Le terme négatif — J R3 \Vi4>f-fo\ 2 dx 
est alors facilement contrôlable en utilisant un nombre fini de fonctionnelles de 
Casimir invariantes. 

- Enfin il reste à étudier la coercitivité du terme 

( / \V x <f) f „ f o\ 2 dx+ / F'{E)\<j> f -f> -V<j)f-p\ 2 dxdv ) . 

Cette dernière implique sans mal une estimation de coercitivité du même terme 
avec Vi à la place de V, pour £ assez grand. Cela revient à étudier la positivité de 
l'opérateur 

A<p=-A x <p+ [ F'{E){(t)-V(t))dv 
Jrz 

agissant uniquement sur le potentiel (j). En remarquant simplement que 

Im(l -P)i_Ker(X>) et Ker(D) = In^D)- 1 , 

on déduit facilement que (0 — V<f) = {h, f } avec h impaire, ce qui permet d'ap- 
pliquer l'inégalité de coercitivité d'Antonov ( 123]) . et de conclure. 
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Ce travail intéressant semble pouvoir se généraliser à des modèles sphériques 
décroissants plus généraux que le modèle de King. La seconde limitation de cette 
méthode, c'est-à-dire le fait de ne considérer que des perturbations à symétrie 
sphérique, semble par contre plus sévère. Un des apports principaux du travail |LMR9j . 
que nous allons maintenant discuter, est de s'être affranchi de cette limitation. 

4.3. Nouvelle approche variationnelle par réarrangement 

Après ce détour par une approche non-variationnelle, nous allons maintenant revenir 
à une approche variationnelle, mais sous un angle nouveau. On voit qu'un défaut de 
l'approche variationnelle par énergie- Casimir est qu'elle semble impuissante à reformuler 
sous forme de problème de minimisation certains modèles stationnaires décroissants. 
Cependant, dans le même temps, l'approche directe par linéarisation semble limitée 
par l'inégalité de coercitivité d'Antonov elle-même et par les difficultés inhérentes aux 
contrôles des termes d'ordre supérieur dans le développement du hamiltonien. 

Le travail |LMR7j constitue une première avancée en introduisant l'idée d'exploiter 
les propriétés d'équimesurabilité du flot : même si la propriété de séparation des états 
fondamentaux n'est pas vérifiée pour le problème de minimisation avec un nombre 
fini de contraintes, l'équimesurabilité de la solution à sa donnée initiale permet de 
prouver dans certains cas une propriété de séparation locale. Finalement dans les travaux 
|LMR8} ILMR9] , Lemou, Méhats et Raphaël résolvent complètement ces contradictions. 
Ils prouvent le théorème suivant dans le cas de modèles sphériques isotropes. 

Théorème 4.2 ( |LMR9j ). - - Soit f° = F(E) > une solution stationnaire continue, 
non nulle, à support compact, du système (JE}-©, pour laquelle il existe E < tel que 
F(E) = pour E > E , F est C l sur ] - oo, E Q [ et F' < sur ] - oo, E [. 

Alors f° est orbitalement stable au sens suivant : pour tous M > et e > il existe 
7] > tel que, pour toute donnée initiale 

Un eS:={g>0, geL 1 ^ L°°(M 6 ), \v\ 2 g G ^(R 6 )} 

telle que 

Wfin-fWmW) <V, Win)<n(f°)+r], ||/i„|U~(R«) < \\f \\ L oo {R e ) + M, 

alors toute solution faible issue de cette donnée initiale vérifie 

Vt>0, / |(1 + \v\ 2 )(f(t,x,v) - f°(x- z(t),v))\ dxdv < s. 

Avant de détailler la preuve, donnons les idées essentielles : 

- en s'inspirant d'idées introduites dans la littérature physique |Oa[ IL-B3[ IWZS[ IAT] 
les auteurs démontrent que le hamiltonien possède une propriété de monotonie par 
rapport aux réarrangements selon l'énergie microscopique ; 

- après ce réarrangement, on est alors ramené à un problème variationnel de mini- 
misation sur le potentiel gravitationnel uniquement, pour lequel la solution station- 
naire est bien un minimum global isolé ; 
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- enfin pour ce problème de minimisation réduit, ils démontrent une inégalité de 
coercitivité d'Antonov généralisée dans ce contexte, et font le lien entre la partie 
radiale de cette inégalité et une inégalité de type Poincaré ; 

- la fin de la preuve est basée sur un argument de compacité pour des suites mini- 
misantes < généralisées » dont le réarrangement selon l'énergie microscopique est 
une suite minimisante pour le problème réduit sur le champ gravitationnel, et la 
compacité est extraite à partir de la coercitivité de l'étape précédente. 

Ce travail met donc à jour une nouvelle structure variationnelle « cachée sous les 
réarrangements selon l'énergie microscopique », pour laquelle l'approche variationnelle 
est bien plus simple et naturelle. Il révèle également le lien entre la coercitivité de 
ce problème de minimisation réduit et un problème d'inégalité fonctionnelle de type 
Poincaré. 

4-3.1. Réarrangement selon l'énergie microscopique. - - Rappelons tout d'abord la no- 
tion classique de réarrangement symétrique (voir par exemple [LL[ Chapitre 3]). Étant 
donné un ensemble A C M 6 mesurable, on définit son réarrangement symétrique A* 
comme étant la boule ouverte centrée en zéro et de même volume que A (pour une 
norme donnée sur M 6 ). Étant donnée une fonction / > intégrable sur M 6 , on définit 
alors son réarrangement symétrique /* comme étant la fonction positive sur R 6 dont 
les ensembles de niveau supérieur sont obtenus par réarrangement symétrique des en- 
sembles de niveau supérieur correspondants de /, ce qui donne la formule suivante par 
intégration par tranche 



r+oo 

f*(x,v)= / l{/>s} ds avec l* A = l A * 
Jo 



La fonction /* est alors radialement symétrique, décroissante, et équimesurable à /. 
Rappelons la propriété élémentaire suivante sur les réarrangements symétriques. 

LEMME 4.3. — Pour f > intégrable sur IR 6 on a 

f*(x, v)\(x,v)\ R 6 dxdv < / f(x,v)\(x,v)\ R 6dxdv 

où \ (x,v)\r6 désigne la norme considérée sur M?. 

Démonstration. — La preuve est très simple, nous la rappelons pour éclairer la suite. 
Pour deux ensembles mesurables A, B C M 6 de volume fini avec \A\ < \B\ (l'autre cas 
étant symétrique), on a 

ï A l B dxdv = \AnB\<\A\ = \A*\ = \A* nB*\= [ l^dxdv. 

Or pour s > donné on a 1*^ v ^ 6<s = 1|(x,v)| R 6<sj e ^ on déduit en utilisant la précédente 
inégalité et en intégrant par tranche 

/ Ai 

Jr 6 Jr 6 
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Puisque J R6 f dx dv = f R6 f* dxdv on en déduit 

/ f*l\( x>v )\ 6>s dxdv < / m6>s àxdv. 

En intégrant finalement selon s G [0, +00 [, on en déduit le résultat. □ 

On introduit maintenant de manière similaire le réarrangement selon l'énergie mi- 
croscopique 

f\v? 
:= (jy + <f>{x, 

d'un potentiel donné </> sur M 3 de la manière suivante. Etant donné un ensemble 4cR 6 
mesurable on définit A*^ son réarrangement selon l'énergie microscopique comme 
étant la « boule d'énergie » ouverte 

A** = {{x,v) | E^x,v)<E A } 

avec Ea choisi tel que \A*^\ = |^4|. Il est facile de voir que 

E^\{(x,v) | E4x,v) < E}\ 

est une bijection de [min£^,0[ sur R + . Étant donnée une fonction / > intégrable 
à support compact sur M 6 , on définit alors f*^ son réarrangement selon l'énergie mi- 
croscopique E^ comme étant la fonction positive sur M 6 dont les ensembles de niveau 
supérieur sont obtenus par réarrangement selon l'énergie microscopique des ensembles 
de niveau supérieur correspondants de / : 



/+00 
l tf>s} dS aVeC ^ = l A 



La fonction est alors une fonction de E^, à support compact, décroissante en l'énergie 
microscopique E^, et équimesurable à /. On a la propriété suivante qui rappelle le lemme 
élémentaire ci-dessus, et dont nous nous servirons par la suite : 

LEMME 4.4. - - Pour f > intégrable à support compact sur IR 6 on a 
f*^(x, v)E ( f ) (x, v) dx dv < / f(x,v)E ( f ) (x,v)dxdv. 

Démonstration. — En raisonnant comme précédemment on a 

1^1,8 dxdi> < / l^fl^dxdw. 

D'où, pour s G [— min £70,0 [, puisque {xv)<s = ^-e^(x,v)<s, en intégrant par tranche 

/ fl E4 ,(x,v)< s dxdv < / f*^l E4 ,{x,v)<sàxd.v 
Je 6 Jmp 

et puisque f R6 f dxdv = f u6 f*^dxdvon en déduit 

/ f^lE^x^ysdxdv < / fl E4 ,( x ,v)>sdxdv. 
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En intégrant finalement selon s G [— min E^, 0[, on obtient 

/ (E^x, v) — min E^) dx dv < f (E^x, v) — min E^) dx dv 

d'où le résultat. □ 

On va par la suite utiliser le réarrangement de / selon l'énergie microscopique créée 
par la fonction f elle-même, que nous noterons / = f*^ f avec, comme précédemment, 
cf)f = — (l/(47r|x|)) * pf. On voit que l'on obtient ainsi une opération de réarrangement 
/ — > f très fortement non-linéaire. Remarquons immédiatement que la solution station- 
naire sphérique est un point fixe de ce réarrangement non-linéaire : f° = (/°)*^/° = f°. 

4-3.2. Monotonie du hamiltonien et hamiltonien réduit. - - On définit la fonctionnelle 

jrW ■.= nu**) + \\\v x <t> - v^||| 2(R3) 

et on montre la propriété de monotonie suivante. 

Proposition 4.5. — Si l'on considère f e S et f — f*^f, alors 

H(f)>J f *(0 f )>H(f) 
avec égalité si et seulement si f = f . 

s 

Ebauche de preuve. — La preuve repose sur le lemme précédent. Par un calcul que nous 
avons déjà fait 

H(f) = U{g) + l -\\V x cf> f - V^JI^a) + f^{^Y + M*)) (/ ~ 9) dxdv 
pour deux fonctions f,g G S, d'où avec g = f : 

W) = Jrttf) + / 6 + M*)) (f - f) ^dv 

et l'on conclut grâce au lemme précédent appliqué à 



L (^2" + 0/(x) ) (f-ÏÏ dxdv ^°- 



Le cas d'égalité se traite en étudiant le cas d'égalité dans l'inégalité de monotonie du 
réarrangement selon l'énergie microscopique. □ 

Il s'avère que la différence Jj* — Jp est facilement contrôlable par des normes sans 
dérivée : 

JAïï - w) > -ii0/iu~( R 3)iir - (/°)iui(R«), 

et l'on ainsi peut se contenter d'étudier la coercitivité de la fonctionnelle Jp. Nous 
noterons J~(4>) '■= J~p{4>) et nous appellerons cette fonctionnelle hamiltonien réduit] 
elle n'agit que sur le potentiel <j). 
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4-3.3. Inégalité de coercitivité d'Antonov généralisée pour le hamiltonien réduit. - 
On va maintenant étudier les propriétés de convexité de J au voisinage de 0/o, et 
ainsi montrer que 0/o est un minimum local de J . On définit un espace de potentiels 
admissibles 

X := <<p E C°(R 3 ) | 0<O, lim0 = O, V x <p E L 2 (R 3 ) , inf (1 + |x|)|0(x)| >0 

[ oo a;6M 3 

On peut alors montrer la proposition suivante. 

Proposition 4.6. — Il existe des constantes cq, 5q > et une application continue 
—> de H l (R 3 ) (l'espace de Sobolev homogène) dans M. 3 telles que pour E X tel 
que 



- 0/0 (• - Z)|| L ^( R 3) + Il V œ - V x (j)f0(- - Z)\\ L 2 



< 5r 



inf 

zGK 3 

alors 

J{4>) ~ J{4>p) > C \\V x (f) - V x 0/o(- - ^)|| i 2 (R 3). 

Ébauche de preuve. — On décompose la démarche en plusieurs étapes. 

Étape 1 : Développement de Taylor. La première phase calculatoire est d'écrire le 
développement de Taylor à l'ordre 2 de J avec un reste contrôlé explicitement : 

J{4>)-J{<Pp) = ^ J D 2 l 7(0/o)(0-0/O,0-0 / o)+O (||0 - 0/o|| Z/ oo (R 3 ) ) ||V :r 0-V :c 0/o||| 2(K 3 ) 

avec le terme d'ordre 1 qui s'annule et 

D 2 J{(f> f0 )(h,h):= [ \V x h\ 2 àx- f \F\E)\[h{x) -{Vh){x,v)] 2 àxàv 

où l'on rappelle que, sans autre précision, E = E^^ = (\v\ 2 /2 + 0/o(a;)), et la projec- 
tion V est définie par 

(\ l 2 

/»3 ( 2 + h° ( x ) - h (y) %) d v 

{Vh)(x,v) :-- 



(\ I 2 \ x l 2 
/r3 (y- + <f>p{x) - 0/o (y) J dy 



C'est un opérateur de projection sur les fonctions de E uniquement. Dans le cas radial 
on retrouve ainsi l'opérateur de projection sur le noyau de l'opérateur v- V x — V x 0/o • V v 
que nous avons déjà rencontré dans le travail [GLj . 

Étape 2 : Inégalité d'Antonov généralisée. Si l'on définit 



£h = -A x h- [ \F'(E)\(1-V)dv 



on obtient 

ï2 



(jCh,h) L 2 m = D 2 J(<f) f o)(h, h) 
et l'on ramène le problème à l'étude de l'opérateur £ ; on a alors la proposition suivante. 
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Proposition 4.7. - - L'opérateur C est positif, c'est une perturbation compacte du 
laplacien sur H 1 (M?) , son noyau est donné par 

Ker(£) = Vect {d xi <f>f>, d X2 (f) f o, d X3 (f) f o } 

et on a donc 

V/i G H\R 3 ), (Ch, h) L2m > c \\V x m 2(R3) - - V ( f hA x (d Xi (f) f0 ) dx] 

c o i=1 \Jr 3 J 

pour une certaine constante co > 0. 

On décompose h en partie radiale et complémentaire orthogonal 

h = ho + h, h G H^ ad (R 3 ), h G (i/ r 1 ad (M 3 )) ± . 

La positivité selon la composante hi est plus simple à traiter car Ilhi = et son étude 
se ramène donc à celle de l'opérateur de Schrodinger A que nous avons déjà étudié plus 
haut dans la preuve de la proposition 13.61 Le noyau Ker(£) de l'énoncé s'en déduit en 
particulier. 

Sur la composante radiale h on a l'inégalité 

y h G HÏ ad (R 3 ), h^O, (Ch, h) L * m > 0. 



Cette propriété peut se démontrer comme dans la preuve de Guo et Lin [GL] que nous 
avons discutée précédemment à la sous-section I4.2[ en montrant qu'elle se réduit à 
l'inégalité de coercitivité d'Antonov démontrée à la proposition 13.51 

Lemou, Méhats et Raphaël proposent une autre preuve intéressante de cette propriété 
et de la proposition 13.5} en faisant le parallèle avec la démonstration d'une inégalité de 
type Poincaré. Ils adaptent la stratégie de preuve de Hormander |Ho-i-ll IHo-i-2j . et 
utilisent une inégalité fonctionnelle de type Hardy. 

Donnons l'idée générale de cet argument. On introduit l'opérateur suivant sur les 
fonctions radiales, exprimé dans les variables E et r : 

Tf(E,r) = 1 -d r f = -4rAf, 

r 2 ^2(E - (fp (r)) r 2 \v\ 

et on vérifie que Vh = implique h = Th pour un certain h. On calcule alors par 
intégration par parties et inégalité de Cauchy-Schwarz (un argument d'approximation 
supplémentaire est nécessaire, que nous n'évoquons pas ici) 

(\ 1/2 
3/ P f o(r) - h \ \F\E)\dxdv 

On montre alors l'inégalité de type Hardy suivante 

/ fp/c(r) + ^M) H \ \F\E)\dxd^\ < I \F>(E)\\Thf dxdv, 
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ce qui, combiné avec l'inégalité précédente, donne 

et conclut la preuve de positivité. L'inégalité de Hardy se démontre en remarquant que 

T 2 h 2 ~, 



{Thf = T (h\h 2 Th 2 \ - ^-h 2 avec h = hji 2 , 



puis 



T2k2 - 3 ^(r) + ^'l pour h 2 = r\2{E-^{r))fl 2 . 



h 2 4r 4 (£-0 f o(r)) 2 

Étape 3 : Traitement du noyau par modulation. On ajuste finalement la fonction de 
translation z$ au moyen d'un théorème des fonctions implicites pour annuler les défauts 
de coercitivité, le., les termes négatifs dans la proposition ci-dessus. □ 

4-3.4- Compacité des suites minimisantes et résolution de la conjecture. - - On peut 
montrer la compacité de certaines suites minimisantes généralisées f n au sens suivant. 



Proposition 4.8. - - Si f n e S vérifie 



sup inf 

n >0 2€K 3 



\\(f) f n - (j) f o(- - Z)IU°°(]R 3 ) + ||Va;0/n - V x (f)fo(- - Z)\\ L 2 



< 5 



et 

lira ||/: - (/ )*||li(r 6 ) = 0, liminf U{f n ) < H(f), 

n— >oo n—ïoo 

alors 

lim / |(1 + \v\ 2 )(f n (x,v) - f°(x - Zé, n ,v))\ dxdv = 0. 

n-too J Re 

Ebauche de preuve. — La preuve est faite en deux étapes. Tout d'abord le contrôle de 
coercitivité précédent implique sans difficultés que 

lim ||V x </>/n - V x (j) f o(- - )\\ L 2i R 3) = 0. 

Ensuite on note f n (x,v) = f n (x + z^ fn ,v) et l'on revient au hamiltonien complet en 
utilisant l'identité 

nfn)-u(f) + hv x 4> f n-v x (f> f 4 L 2 m = f E^jû-f) dxdv. 

A partir des hypothèses et de la convergence déjà démontrée, on a 

lim sup (H(f n )-H(f°)) = lim sup (H(f n )-H(f )) < 0, 



\V x (j)fn ~ V x 4>f0\\ L 2 {R3) = \\V x (j)fn ~ W x (j)f0{- ~ ^ /n )||L2( R 3) -> 0, 



et on en déduit 



lim sup / (/„ - f) dx dv < 0. 
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L'hypothèse (f n )* — > (/°)* implique par ailleurs 

lim / ^ (/°-^)drdT; = 



d'où l'on déduit 



lim sup / E^ (f n - f^ f ° ) dx dv < 0. 

n— S-oo 



De par la monotonie du réarrangement j RIS E^ i0 {f n — f n ^ f0 )dxdv > 0, on en déduit 
finalement 

lim / E hQ {fn-fn i0 )àxàv = Q. 

Il suffit ensuite de montrer que la saturation de cette inégalité de réarrangement, com- 
binée à l'hypothèse (f n )* —> (f )* implique que 

lim ||/n-/°|Ui(R«) = 0. 

En combinant cette convergence avec l'hypothèse de limite supérieure sur le hamiltonien 
ainsi que la convergence du potentiel, on obtient finalement la convergence de l'énergie 
cinétique 

lim / \v\ 2 f n dxdv= / \v\ 2 f°dxdv 
ce qui conclut la preuve de la proposition 14.81 □ 



La fin de la preuve du théorème 14.21 se fait ensuite en combinant la proposi- 
tion 14.81 l'inégalité d'interpolation de la proposition 12.14 ainsi que la contractivité du 
réarrangement symétrique ||/* - (/°)*||z,i(K6) < 11/ - 



5. CONCLUSION ET PROBLÈMES OUVERTS 

Ce problème de stabilité des galaxies est un exemple intéressant de recherche 
mathématique nourrie par une question concrète posée par la physique théorique. Nous 
essayons pour terminer de soulever quelques questions ouvertes d'ordre mathématique, 
en suggérant des liens avec d'autres travaux. 

Tout d'abord, la première question naturelle du point de vue de la pertinence physique 
des résultats est de quantifier la taille du voisinage de stabilité orbitale. Cela paraît 
maintenant une tâche plus abordable avec la nouvelle théorie de Lemou, Méhats et 
Raphaël ; il s'agit essentiellement de rendre explicites, ou tout au moins constructives, 
les constantes de coercitivité dans les inégalités fonctionnelles utilisées. 

L'autre question naturelle est de sortir du cadre strictement monotone pour la solu- 
tion stationnaire f°(E). Par exemple, nous pouvons déjà nous demander si, au niveau 
des solutions stationnaires, localement au voisinage d'une solution orbitalement stable, 
il est possible de démontrer un théorème de paramétrisation bijective des solutions sta- 
tionnaires par les conservations du système, dans le même esprit que le travail récent 
de Choffrut et Sverâk |CSj sur l'équation d'Euler incompressible en dimension 2. 
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Cependant, nous pourrions nous attendre plus généralement, au niveau dynamique, 
à la stabilité orbitale autour d'une solution stationnaire « presque » monotone, et donc 
proche des solutions orbitalement stables que nous avons étudiées. Une première tâche 
serait ici de clarifier au niveau mathématique les instabilités créées par des perturbations 
non radiales de modèles sphériques anisotropes. 

Les méthodes variationnelles semblent néanmoins trouver leur limite, et cela soulève 
la question de revenir à nouveau à une approche directe par linéarisation. Cela nous 
amène également à faire une autre remarque importante sur les travaux que nous avons 
présentés : ceux-ci n'utilisent pas la dynamique proprement dite de l'équation, mais 
uniquement ses invariants. Même si les résultats obtenus sont dynamiques, le cœur 
conceptuel de ces méthodes n'utilise pas la dynamique. C'est une force de ces approches, 
qui leur confère une grande robustesse pour traiter des modèles généraux et manipuler 
des solutions très faibles, mais c'est également une faiblesse dès lors que l'on sort d'un 
cadre parfaitement variationnel. 

Par conséquent, il serait intéressant d'explorer les dialogues possibles avec les résultats 
de stabilité non-linéaire obtenus dans |MV] . Ces derniers résultats utilisent des espaces 
fonctionnels très réguliers, c'est donc en ce sens l'extrême opposé des méthodes de 
stabilité orbitale que nous avons présentées. En particulier, cela implique qu'il faut 
« traquer » les oscillations du système dans les estimations de régularité, alors que les 
espaces de Lebesgue utilisés dans les théories de stabilité orbitale ne « voient » pas les 
oscillations de la variable de vitesse produites par le mélange de phase. Les hypothèses 
sur les données initiales sont donc bien plus fortes, mais cela permet également d'obte- 
nir une information plus précise sur le comportement asymptotique du système, et de 
montrer la stabilité non-linéaire autour de solutions stationnaires linéairement stables, 
mais qui ne vérifient pas un problème variationnel. 

Les résultats de |MV] s'appliquent à l'équation de Vlasov-Poisson gravitationnelle, 
mais uniquement dans le cas non physique d'un domaine périodique en espace. Le cas 
d'un système auto-gravitant qui « crée sa propre géométrie » et son propre confinement 
au cours du temps représente un défi probablement difficile mais aussi très intéressant 
pour cette approche par linéarisation : à l'inverse de |MVj où « l'amortissement Lan- 
dau » produit une convergence vers zéro du champ moyen asymptotiquement, on ne 
connaît plus ici à l'avance quelle doit être la limite du champ moyen lorsque t — > +00. 
Un défi conceptuel similaire se pose pour l'équation d'Euler incompressible en dimen- 
sion 2 ainsi que pour l'équation de Vlasov-Poisson avec un confinement magnétique. 

Pour terminer, nous mentionnerons le problème intéressant, mais probablement pour 
le moment hors d'atteinte tant que les questions précédentes et la limite de champ moyen 
ne sont pas mieux comprises, de faire le lien entre les résultats de stabilité obtenus pour 
l'équation de Vlasov-Poisson gravitationnelle, et les résultats de stabilité asymptotique 
(on pense par exemple à la théorie KAM) pour le problème à iV corps, dans la limite 
N -»■ 00. 
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